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Band structure properties of 182 Os are investigated through a particle number and angular momen- 
tum constrained generator coordinate(GCM) calculation based on self-consistent three-dimensional 
cranking solutions. From the analysis of the wave function of the lowest GCM solution, we confirm 
that this nucleus shows a tilted rotational motion in its yrast states, at least with the present set 
of force parameters of the pairing-plus-quadrupole interaction Hamiltonian. A close examination of 
behavior of other GCM solutions reveals a sign of a possible occurrence of multi-band crossing in the 
nucleus. Furthermore, in the course of calculations, we have also found a new potential curve along 
the prime meridian on the globe of the J = 18h sphere. Along this new solution the characters of 
proton and neutron gap parameters get interchanged. Namely, A p almost vanishes while A n grows 
to a finite value close to the one corresponding to the principal axis rotation(PAR). A state in the 
new solution curve at the PAR point turns out to have almost the same characteristic features of an 
yrare s-band state which gets located just above the <?-band in our calculation. This fact suggests 
a new type of seesaw vibrational mode of the proton and the neutron pairing, occurring through a 
wobbling motion. The mode is considered to bridge the g-band states and the s-band states in the 
backbending region. 
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I. INTRODUCTION 

A limitation of a framework of the nuclear mean field 
approach whose validity is supported by the success of 
the idea of an independent particle motion has become 
one of the central matter of concern in recent years in 
nuclear physics. When we describe, for instance, high- if 
states in a 7-soft nucleus or a large amplitude collec- 
tive motion(LACM), such as a fission or a process which 
leads to break up into fragmentation, which has recently 
attracted much attention through the molecular dynam- 
ics analysis, an assumption of stiffness against deforma- 
tion is no longer valid. Furthermore, in the nuclear halo 
or in a skin structure only a small number of nucleons 
come into play in the region, and therefore the notion 
of a formation of a stable mean field will hardly be ac- 
ceptable in such exotic situations. For the description 
of these new topics one thinks of the two representative 
theories: the time dependent Hartree-Fock(TDHF) and 
the generator coordinate method Q. They demonstrate 
their full ability in description of such situations. Partic- 
ularly, the GCM is fully quantum mechanical and it can 
be applied to wide range of problems. Especially in the 
situations where several potential minima come into play 
in the motion. For instance, LACM; in the situations of 
including a super deformation or a hyper deformation, it 
will be one of the promising methods because it expresses 
states as a superposition of many different product wave 



functions and can describe tunneling effects among the 
potential minima. The basic idea of this method which 
goes beyond the conventional mean field theory, may well 
continue as one of the reliable frameworks in the coming 
century. 

Since the nuclei in the region of osmium and tungsten 
generally exhibit softness against the 7-deformation Q 
they provide a nice testing ground for an approach which 
goes beyond the mean field theory. We have been carry- 
ing out a systematic calculation based on a self-consistent 
three-dimensional cranked Hartree-Fock-Bogoliubov(3D- 
CHFB) theory in 182 0s and have found a solution which 
suggests the existence of a tilted axis rotating(TAR) state 
Q . The resulting potential curve along the prime merid- 
ian on the globe of J = 18h sphere has two TAR minima 
which locate symmetrically with respect to an equatorial 
plane. Thus the GCM calculation by taking into account 
the north and the south latitude as a generator coor- 
dinate to describe excited states is a highly interesting 
problem. 

Our previous attempts based on the GCM calculations, 
however, suffer from some difficulties [0. The first point 
was that in our calculation somewhat large energy split- 
tings among the GCM solutions were obtained. Since we 
wish that our GCM solutions should correspond to the 
observed level spacings above the yrast states, the calcu- 
lated energy spacings among these solutions should stay 
within the range of 10keV~100keV. The second difficulty 
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in our calculation was that we hardly obtained any stable 
GCM solutions in the high angular momentum region. 

One of the problems in our approach is that we lose 
much information concerning the high frequency part of 
wave functions when we solve the GCM equation. Be- 
cause in solving the GCM equation we have discarded 
the eigenstates of the norm overlap kernel whose eigen- 
values are negligibly small. These states bring serious 
numerical errors when the equation is numerically solved 
p2[ . Therefore, they are considered to be unfavorable. 
The other possible problem is that although the TAR 
potential spread non-locally over a surface of a globe 
with a constant angular momentum, we restricted the 
Hill- Wheeler integration to one-dimensional coordinate, 
that is we restricted the calculation only along the prime 
meridian. These drawbacks of our approach will be reme- 
died by introducing the the angular momentum projec- 
tion techniques [|| and the two-dimensional Hill- Wheeler 
integration |^] . Since the programs for carrying out such 
calculations will become a project of inevitably exhaus- 
tive investigation, more basic study about the character- 
istic features of GCM solutions at this moment must be 
helpful. 

Generally, in the process of solving the GCM equation 
the expectation values of the angular momentum and the 
particle number may get shifted from the values of the 
CHFB constraints. For instance, the amplitudes of lower 
angular momentum components present in the cranking 
wave function will increase their probabilities in such a 
way to reduce the total energy of the system. To keep the 
correct configuration of wave functions we restrict the ex- 
pectation values of the angular momentum and particle 
number to their given values during the process of solv- 
ing the GCM equation. Applying the constrained GCM 
method @ to the TAR potential curve to obtain the im- 
proved solutions is one of the main motivations of the 
present work. 

By solving the constrained GCM equations in case of 
182 Os at J = I8T1 we have found that the wave func- 
tions have symmetric or anti-symmetric properties with 
respect to the generator coordinate. This means that the 
broken symmetry of signature in TAR wave function has 
been recovered through the GCM procedure. Analyzing 
the wave function we have found that this nucleus must 
be a tilted-rotor in the present parameter set. Further- 
more, we have found the existence of the PAR compo- 
nents in the excited GCM solutions. This result is very 
interesting because we have solved the GCM equations 
based on the TAR potential which shows the minimum 
points located at a distance from the principal axis. This 
fact exhibits a sign of existence of a different band struc- 
ture in the GCM solutions at J = 18h state. We will 
discuss a possibility of multi-band crossing in this nu- 
cleus. Some preliminary results of our study are already 
presented in a short article M. 

In the next section 2 a brief sketch of the formalism 



is presented. Results and discussions are presented in 
section 3. Finally section 4 contains a brief summary 
and conclusions of the present investigations. 



II. CALCULATIONAL METHOD 

In this section we describe briefly the framework of our 
model. Details of the calculational methods for the 3D- 
CHFB and for the Hill- Wheeler equation are presented 
in our previous works p|,|9|,|l0[ . 

The model Hamiltonian used in the present study is 
the pairing plus Q-Q force given by, 
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where represents the single particle energies and cj 
and Ci are the creation and the annihilation operators 
for particle states | i) . Conditions Cj | 0) = hold for all 
i, and | 0) stands for the true vacuum. The parameters k 
and g T are the corresponding interaction strengths. The 
quadrupole operator is defined, 
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This Hamiltonian is very simple for calculation but still 
effective in describing essential features of the nuclear 
deformation and the pairing correlation. 

In the present work the Hill- Wheeler equation is mod- 
ified to include the constraint on the expectation values 
of the angular momentum and the particle number op- 
erators. As usual, first the norm overlap kernel is diago- 
nalized, 
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where n(i/j,tl>') is Hermitian matrix defined as, 

n(^V') = W)W)>- (5) 

The wave function Q(ip) is the solution of the following 
3D-CHFB equation with eight constraints, 
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where fx k , £,k and A r are the Lagrange multipliers which 
are adjusted to satisfy the following constraints, 
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and 
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where B k 's express three off-diagonal elements of the 
mass-quadrupole tensor in Cartesian-coordinate repre- 
sentation, 
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We know from our previous study that a calculation 
along the prime meridian always gives lower energy when 
compared with the calculation performed on the whole 
surface of the globe ||. Therefore, in the present study 
we restrict the generator coordinate only to the latitude 
if) along the prime meridian. Thereby the CHFB equa- 
tions become now only two-dimensional. The angle tp is 
measured with respect to the cc-axis so that Ji = J cosip, 
Ji = and J 3 = J simp. This also helps us not to face 
any problem in choosing the phase in the square root 
which will appear in the calculation of eq.(5) pT| . By 
utilizing the eigenfunctions u k (ip) of the overlap kernel, 
eq.(4), we can construct the orthogonal sets which are 
called the natural states |l^ ], 
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Generally, in the process of solving the GCM equa- 
tion the expectation values of angular momentum and 
the particle numbers will somewhat get shifted from the 
values of the CHFB constraints. To keep the appropriate 
wave functions we restrict the expectation values of the 
scalar of squared angular momentum vector J 2 and par- 
ticle numbers to their correct values during the process 
of solving the GCM equations. Therefore we solve the 
following constrained Hill- Wheeler equation, 
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The Lagrange multipliers Af* and in eq.(12) con- 
trol the expectation values of particle numbers and the 
square of angular momentum in the GCM wave functions 
and are different from those appearing in eq. (6) . The su- 
perscript "a" indicates the index of the eigenstates. The 
matrices H k ^ k *, N k k , and j\ k , are given as follows, 
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Detailed expressions for the matrix elements appearing 
in the integration are given in Ref. [f|. The eigenstates 
of our GCM equation are given by, 
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where g k °^ are obtained as solutions of eq.(12). From 
eq.(ll) we find the corresponding weight function, 
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Then, finally the physical quantities are calculated using 
the matrix elements in eq. (13-15) and the wave function 
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as follows, 
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III. NUMERICAL RESULTS AND DISCUSSIONS 

A. Solution of 3D-CHFB equations 

Recently a systematic calculation based on the 3D- 
CHFB theory revealed the possible existence of a tilted 
axis rotating potential in 182 Os nucleus ||[|. In our pre- 
vious calculations, however, the force parameters were 
chosen so as to reproduce just the overall features of low- 
angular momentum states of nuclei in this mass region. 
The observed first backbending in A~180 region is ex- 
pected to be due to «i3/2 neutron alignments, whereas in 
our earlier calculation protons were playing the role. The 
reason was attributed to our slightly too strong Q-Q force 
parameter and also to rather strong pairing interaction 
for neutrons. In the present calculation we reduced the 
Q-Q force parameter and the strength of neutron pairing 
interaction so as to make rotation alignment of neutron 
occur in the PAR states, which is considered to be yrast 
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in 8 Os. Our new force parameters are chosen so as to 
give as (3 = 0.280, A p = 1.040MeV and A„ = 0.925MeV 
in the ground state. Moreover, we have attained con- 
siderable high precisional calculation this time. We have 
increased the number of iteration to achieve better con- 
vergence when solving by the steepest descent method. 
We have fixed an increment of angular momentum as 
A J = Q.lh in the calculation along the yrast line up to 
J = 3QTi and set the maximum iteration for convergence 
as 300 times and the minimum iteration as 50 times for 
each step. In our previous calculation, those iteration 
times were taken as 120 and 30, respectively. The size of 
the model space and the single particle energies q in the 
Hamiltonian are not changed from our previous study. 
They are shown in Table 1 together with quantum num- 
bers such as isospin, radial node number, orbital angular 
momentum and total angular momentum [fl3"[ . 

In a similar manner as the previous papers, we ob- 
tained PAR solutions by cranking up along x-axis per- 
pendicular to the symmetry axis of the non-rotating state 
and also by cranking down along the axis. The two bands 
cross at J ~ 20fr in the present case. At the crossing 
point in the g-band the size of proton gap parameter is 
A p = 0.87MeV and the neutron gap parameter almost 
vanishes. The characteristic features of the yrast line are 
indicated in Fig.l. 

Model space and single particle energies 



36 valence particles and 70 core particles are considered 
for the neutron. 



proton 


orbits 


energy (MeV) 


I/7/2 


8.06435 


0«13/2 


7.55675 


O/19/2 


7.23480 


2si/ 2 


3.87120 


14/2 


3.59468 


Oftll/2 


2.16538 


14/2 


1.29040 


057/2 


0.64520 


039/2 


-3.50251 


neutron 


orbits 


energy (MeV) 


0«H/2 


7.55740 


109/2 


6.74629 


2P1/2 


4.40579 


1/5/2 


3.39191 


2R3/2 


3.02322 


0«13/2 


1.56626 


Ohg/2 


0.82954 


1/7/2 


0.16591 


Ohn/2 


-4.23988 



Table 1: Model space and the single particle energies in 
the present calculation. In the space 36 valence particles 
and 40 core particles are taken for the proton and also 
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FIG. 1. Fig. la shows the behavior of energy curves of the 
g-band and the s-band as well as the experimental value. The 
two bands cross at J ~ 20ft in our calculation. The energy 
value 211keV indicated in the figure shows the energy differ- 
ence between the J = 187i yrare state and the yrast state. 
Fig. lb shows the behavior of the gap parameters for each 
band. The change of 7-deformation is indicated in Fig.lc. 

Fig. la shows the behavior of energy curves of the g- 
band and the s-band as well as the experimental value 
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taken from Ref. 0. In our previous calculation M both 
bands cross at J = 14fi, due to the proton alignment. In 
the present case, after the band crossing at J = 20h the 
(/-band continues up to J = 25h state. We can not obtain 
any <?-band components above this state and after the 
intensive iteration processes the CHFB solution finally 
converges to a next state which belongs to the s-band. 
In the same manner, the s-band obtained by cranking 
down calculation continues down to J = 14h state. When 
decreasing the angular momentum by a unit O.lh from 
this state the solution converges to a new state in the g- 
band after many iteration processes. In Fig. lb we show 
the behavior of the gap parameters. The strength of the 
neutron gap parameter both in the g-band and the s- 
band are much reduced when they are compared with our 
previous calculation [0]. In the present calculation, the 
proton gap parameter almost keeps its value until the end 
of the <?-band. The change of 7-deformation is indicated 
in Fig.lc. The 7-deformation changes drastically at the 
band crossing point. The value in the g-band is 7 S = 
7.36° and in the s-band is 7,5 = —10.19° at the point. 

We tilted rotating axis along the prime meridian, which 
intersects with z- and :r-axis on the sphere of J — 18?i. 
As in our previous calculation, here too the existence of a 
stable TAR state is confirmed in the region of ip = 20° ~ 
30° as shown in Fig. 2a. 

The TAR minimum takes place at tp = 24°, and 
its depth is Vtar = — 0.166MeV. We can find an oc- 
currence of phase transition between the PAR and the 
TAR point at about tp = 10°. Starting from the PAR 
the potential curve ends at this point with its energy 
Vfi rs t = 0.189MeV. After this point, a new potential 
curve with different configuration appears which contin- 
ues up to tp = 42°. In the region between tp = 10° and 
tp — 42°, in which the TAR appears, the proton gap pa- 
rameter almost vanishes and the neutron gap parameter 
recovers its original strength of the ground state value as 
shown in Fig. 2b. This fact implies that the TAR is gen- 
erated due to a collapse of the proton gap parameter in 
our calculation. Thus the phase transition is considered 
to be caused by a crossing of a band with A„ = and 
A p to an another band with A„ 7^ and A p = 0. 

We encounter the next phase transition point at 43° 
with V secon d = 0.470MeV. At this point the configura- 
tion of wave function also changes drastically. That is, 
the proton gap parameter recovers its strength and op- 
positely the neutron gap parameter reduces its value to 
about 70% of the value at the TAR region as shown in 
Fig. 2b. In the region between tp — 44° and tp = 52° 
both the proton and the neutron gap parameters take 
about A p ~ A„ ~ 0.5MeV. Calculation above this re- 
gion, where we did not explore in our previous GCM 
study, the potential curve shows relatively smooth behav- 
ior and it goes down into negative value at the north pole. 
The energy value at the point is V n - po i e = — 0.167MeV. 
This means that the states around the north pole are 



lower in energy compared with the PAR state. The wave 
function of this state consists of the positive parity H3/2 
neutron which align their angular momenta to the sym- 
metry axis and contribute to K = 8 band state. At this 
moment, we do not have any experimental evidences for 
the TAR minimum or the existence of energetically lower 
states near the north pole. 
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FIG. 2. Fig. 2a shows the energy curves along the prime 
meridian on the globe of J = 18h sphere of 182 Os. The solid 
curves (indicated "forward") are obtained by increasing the 
north latitude angle starting from tp = 0° up to tp = 90° . The 
two phase transitions are observed at tp = 10° and tp = 42°. 
The two dashed curves show the another 3D-CHFB solutions 
(indicated "backward" ) which are obtained by decreasing the 
angle. Along this new solutions the characters of the proton 
and the neutron gap parameters mutually get interchanged 
as seen in Fig. 2b. All the energy curves are mirror symmetric 
with respect to an equatorial plane of the globe. Fig. 2b shows 
behavior of the pairing gap parameters for the proton and the 
neutron along the curves indicated in Fig. 2a. We can see that 
the TAR potential is caused by a sudden decrease of proton 
gap parameter and the simultaneous sudden restoration of the 
neutron gap parameter. Near the PAR point the proton gap 
parameter A p almost vanishes and the neutron gap parameter 
A n recovers its value. 

The GCM study based on this full potential curve will 
certainly give some clues to understand the new dynamics 
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characteristic in this mass region. Unfortunately, our 
previous calculation were restricted to a region up to the 
second phase transition point 0. 

On the other hand, we show in Fig. 2a the existence of 
another 3D-CHFB solutions along the prime meridian in 
this nucleus [H. The first solutions have been obtained 
by a calculation starting from the point ip = 30° back to 
the point ip = 0°. The next solutions starting from the 
point ijj = 62° back to the point ip = 32°. Along this 
new first solution the characters of proton and neutron 
gap parameters mutually get interchanged as shown in 
Fig. 2b. Namely, A p almost vanishes while A n grows up 
to a finite value at the PAR point. Therefore this new 
solution has an s p -band character. At the PAR point, 
the excitation energy and the gap parameters A p and 
A„ of this new solution takes the close value as that of 
the yrare s-band state at J = 18h. The excitation en- 
ergy of the new solution at the PAR point is E t ut-back ~ 
164keV whereas the energy difference between an yrast 
g-band state and an yrare s-band state at J = 18h is 
Eg^ s ~ 211keV in our calculation(see Fig. la). The en- 
ergy difference between the states in the excitation en- 
ergy is about 47keV. At this moment, we think that this 
difference arises due to still insufficient convergence dur- 
ing our iteration process. If we start a tilting calculation 
along the prime meridian from the E — 211keV J = 187i 
yrare state, a resulting potential curve just after the first 
phase transition point at ip — 10° (see Fig. 2a) appears 
very close to the original potential curve which is ob- 
tained by starting from the yrast J = 187i PAR state. 
The new state at the PAR point can naturally be con- 
sidered as a member of an yrare s-band states, but more 
precise study will be necessary about this point. The 
value of gap parameters at the PAR point in the new so- 
lution are almost the same as that of the corresponding 
yrare state in the s-band. 

In the present calculation the s-band solutions in the 
PAR calculation continue down to a state with angular 
momentum J = 14?i. Along the solution curve starting 
inversely from the J = 30?i the neutron gap parameter 
gradually increase up to A„ = 0.76MeV at J = 16h. The 
7-deformation changes drastically between the yrast and 
the yrare states at J = 18?i. The value of the yrast state 
is 7 = 5.88° and that of the yrare state is 7 = —9.77°. 
From these results, we can say that a new type of seesaw 
vibrational mode of the proton and the neutron pair- 
ing, occurring through wobbling motion accompanying 
the large 7-oscillation, can bridge the yrast state in the 
g-band and the yrare state in the s-band in the back- 
bending region of 182 Os. 

B. Solving GCM equations 

We have solved the constrained Hill- Wheeler equation 
eq.(12) for 182 0s at J — 18ft. in which the north and the 



south latitude angles are taken as the generator coordi- 
nate. The integration has been carried out in the range 
of ip = —90° ~ 90°. In numerical calculations we replace 
the integration over the range into a discrete summation 
with a step of 4° starting from i/j — —88° to 88° in which 
totally 45 integration points are contained. The size of 
the step should be chosen carefully. In the GCM calcula- 
tion, we first diagonalize the overlap kernel to construct 
a complete orthogonal set of the GCM wave function. If 
we set the size to smaller than 4° in the calculation, some 
eigenvalues of the norm overlap kernel turn negative be- 
cause of an over-completeness of the overlap kernel. In 
our calculation all eigenvalues of the norm kernel stay 
positive and the smallest one is in the order of ~ 10~ 7 
which stays in the range of numerical error. 

In this study we have carried out two types of calcula- 
tions. The first type of calculation is the numerical solu- 
tions without any constraints whose solutions are shown 
in Fig. 3a. The second one is the solutions in which the ex- 
pectation values of the squared angular momentum and 
the particle number operators are constrained. The re- 
sults are shown in Fig. 3b. The lowest ten excited states 
are shown in both the figures. The abscissa of the figures 
are the number of eigenstates taken into the calculation 
of eq.(12) which are obtained from a diagonalization of 
the norm overlap kernel. The eigenstates are taken in the 
decreasing order of the size of eigenvalues of eq.(4). 

In Fig. 3a, we can see that the solutions do not show any 
stable plateau in the whole range of abscissa coordinates. 
The situation is not improved at all when compared to 
our previous works. We only find relatively small plateau 
in the region of lower excitation energies within a range 
of abscissa coordinate 35 ~ 40. The lowest two solutions 
in the region are almost degenerate and the next higher 
state lies more than 700keV higher in the energy position. 
In the previous calculation we also obtained relatively 
large energy splittings among the GCM solutions and we 
attributed the reason to relatively deep potential of the 
TAR. The value was V^ R = -0.367MeV whereas in the 
present case the value is Vtar = — 0.166MeV. The depth 
is much smaller this time and the value is less than half 
of the previous one. From this fact, we can consider that 
only the shallowness of the TAR potential may not be 
enough to suppress the size of energy splittings among 
the states. 

In Fig. 3b, we can find an existence of relatively large 
and stable plateaus in the whole range of the abscissa co- 
ordinate. The energy splittings among the states are in 
the range of lOkeV ~ lOOkeV. These values are reason- 
able for an interpretation of any excited states observed 
in this mass region. At least the lowest seven states con- 
sists relatively stable plateau within the range of abscissa 
coordinate 30 ~ 35 (here after we refer this region as regu- 
lar) . From the present investigation, we understand that 
the conservation of the square of angular momentum and 
the particle number can very much increases the stability 
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of the GCM solutions. Usually a large angular momen- 
tum fluctuation is included in the wave function resulting 



from a cranking model. Since this value is A J 2 
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in our case, we restrict the expectation value of J to 
<~ 420h 2 for studying the J = 18h state. In solving the 
Hill- Wheeler equation the expectation value of the to- 
tal angular momentum does not remain constant when 
increasing the number of states but drifts to lower an- 
gular momentum value to reduce the total energy of the 
system. This situation can be clearly seen in Fig. 4a. 
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FIG. 3. Several lowest solutions of Hill- Wheeler equations. 
In the figures the abscissa indicates the number of solutions of 
eq.(4) considered in solving the equations. Fig. 3a shows the 
solutions without any constraints and Fig. 3b indicates the 
results with particle numbers and the square of angular mo- 
mentum are constrained (see eq.(12)). Since the wave function 
in the cranking model suffers from a large angular momentum 
fluctuation, we took the constrained value of the squared an- 
gular momentum value as J 2 ~ 4207i 2 in the calculation for 
the J = 18h state. We can see from the figure that the con- 
straints very much contribute to stabilize the GCM solutions. 
In Fig. 3b we indicate the point at which GCM weight func- 
tions eq.(17) are calculated(see Fig. 6). 

The solutions without any constrains do not keep the ex- 
pectation value of angular momentum. Contamination of 
lower angular momentum components into our GCM so- 
lutions may force the excitation energy of each solution 
to lower values. This is the main reason for obtaining 



the unstable GCM solution when solving with no con- 
straints (Fig. 3a). The angular momentum constraint will 
not be necessary if good angular momentum states are 
used for the GCM calculation. 
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FIG. 4. Comparison of behaviors of the constrained val- 
ues when solving the Hill- Wheeler equation with and without 
the constrained terms. Fig. 4a indicates the difference of J 2 
between the two results. We can see that the decrease of 
the expectation value of angular momentum in Fig.4a corre- 
sponds to the instability of the Hill- Wheeler equation shown 
in Fig. 3a. Fluctuation for the neutron numbcrs(Fig.4b) and 
for the proton numbers(Fig.4c) are small in the present study. 

In our present study constraints for the particle num- 
bers do not appear to be essential. Following the above 
discussions, contamination of more particle number com- 
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ponents in the wave function may lower the total energy 
of the system. But this tendency may be suppressed by 
the contamination of lower angular momentum compo- 
nent in the wave function. In our calculation, the fluctu- 
ations of particle numbers in the solutions of GCM equa- 
tions without the constraints are negligibly small, that is, 
the fluctuation for the neutron is less than 0.56%(Fig.4b) 
and for the proton the value is almost unchanged(Fig.4c). 

Fig. 5 shows the probability distribution of the GCM 
solutions. In the plane the radius of each circle indicates 
the size of eigenvalue of the equation of norm overlap 
kernel, eq.(4), which express the probability of each so- 
lution. These solutions are taken into account in solving 
the GCM equations. For instance, the slope of a line 
for an energy versus J 2 plot (Fig. 5a) indicates a value 
for the chemical potential M « in the eq.(12). In our 
calculation this value becomes 6.20 x 10~ 3 within the 
regular region. From the figure we can see that the 
solutions are distributed approximately along this line. 
The most probable GCM states localize around the re- 
gion J 2 = 400?i 2 <~ 450H 2 and their energies are about 
E' 42MeV. 

Fig. 5b and Fig. 5c show the distribution for the GCM 
solutions in the energy versus particle number plane. 
We can see the most probable states locate around the 
value of 36 both for the protons as well as for the neu- 
trons. The chemical potentials in the regular region are 
A^ 1} = 4.35MeV for the proton and a! 1} = 3.66McV for 
the neutron. These values are almost constant within the 
region. 



C. Characteristic features of GCM solutions 

In Fig. 6 we show behavior of the weight functions ex- 
pressed by eq.(17) for several of the GCM solutions which 
are obtained from the constrained Hill- Wheeler equation 
(12). In the calculation, the number of states obtained 
by eq.(4) which are taken into account in solving eq.(12) 
are "32" . That point locates within the regular region 
and is indicated in Fig. 3b. From Fig. 6 we confirm that 
the GCM solutions have symmetric or anti-symmetric 
properties with respect to the north and the south lati- 
tudes. This means that the broken symmetry of signature 
in TAR wave function has been recovered through the 
wobbling motion. All the GCM solutions split into two 
groups, one is a symmetric group and the other is an anti- 
symmetric group. In the present case, the first, fourth, 
fifth, sixth and eighth solutions belong to the symmetric 
group and the second, third and seventh solutions belong 
to the anti-symmetric group. 
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FIG. 5. The probability distribution of the GCM solutions. 
In the plane the radius of each circle indicates the size of 
eigenvalue of the equation of norm overlap kernel. The size 
express the probability of each solution. From the slope of 
this plot we can reduce the chemical potentials fx^, A^, 1 ' and 
^ which appear in eq.(12). 
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FIG. 6. The weight functions for the several lowest solu- 
tions of constrained Hill- Wheeler equation (eq.(12)) at the 
point indicated in Fig. 3b. We can see that the GCM solutions 
have symmetric or anti-symmetric properties with respect to 
the north and the south latitudes. This means that the bro- 
ken symmetry of signature in TAR wave function has been 
recovered through the wobbling motion. 

The lowest solution expresses that this nucleus in the 
yrast J = 187i state is a tilted- rotor at least for the 
present set of parameters. The tilt-angle is in the range 
of 12° - 22°. The most probable angle is ip = 12°. 
In Fig. 2a we show the TAR potential curve which has 
a minimum point at ip — 24°. There is no guarantee, 
however, that a nuclear rotational motion will take place 
about an axis indicated by this angle. In the present ap- 
proach, only the GCM solutions based on the potential 
curve can properly describe the characteristic features of 
nuclear collective motion, such as a tilted axis rotation. 
This point was not clearly discussed in our previous pa- 
pers |^|,[)|. Our lowest solution is symmetric with re- 
spect to ip = 0° and is thus characterized as gerade state. 
This solution also exhibits a vestige of a PAR character 
since it shows a small amplitude at ip = 0°. As already 
discussed in section 3.1, the PAR state at J = 18h in this 
nucleus is characterized by the almost vanishing neutron 
gap parameter and the still surviving proton gap parame- 
ter. Therefore we expect that this small PAR component 
of the state shares a neutron aligned s-band character. 
We call such state as an s„-band state. These PAR char- 
acteristic features become remarkable in other solutions 
obtained at higher excitation energies as discussed below. 

The second and the third solutions also show the fea- 
ture of the tilt states. The most probable tilt-angle for 
both states is also ip = 12°. The second one, however, 
shows the tilt characteristic feature also at ip = 8°. These 
two solutions have no amplitudes at the origin and show 
anti-symmetric character with respect to ip = 0°. Thus 
they are grouped into the ungerade states. 

The fourth and sixth states are the gerade states and 
indicate a nature of the PAR state because they show cer- 
tain amplitudes at ip — 0° . This result is very interesting 
because we have solved the GCM equations based on the 



TAR potential curve which has the minimum points lo- 
cated at a distance from the principal axis. Since the 
state with an angular momentum J = 18h is located in 
the neighborhood of the band crossing point in 182 0s, we 
can naturally expect states with an s-band character in 
the GCM solutions. If above symmetric solutions are the 
member of an s„-band, understanding of the backbend- 
ing phenomena in this nucleus becomes quite interesting. 
P.M.Walker et al. suggested that the backbending 
is expected to occur due to crossing of a <?-band with 
a f-band and an s-band is not expected to come into 
play in their study. However, if such a PAR-like solution 
will come into the lowest state, that is an yrast, in the 
higher angular momentum region, the situation warrants 
a somewhat different scenario. The possibility of s-band 
crossing after the i-band crossing emerges as one of the 
backbending mechanisms in this nucleus. Therefore, we 
have found a sign of multi-band crossing in this nucleus. 

The fifth solution also turns out to be a tilt state whose 
tilt angle is tp = 20°. Since the amplitude of the weight 
function between these two tilt points is almost negligi- 
ble, this state may be considered as almost a pure tilted- 
rotor. 

By looking at the higher angle region in Fig. 6, we can 
find that all the lowest seven solutions show certain am- 
plitude in the range of ip = 44° ~ 66°. This result sug- 
gests an existence of a high- if states (K = 12 ~ 16). 
Further interesting features in the results are that the 
fourth and the fifth solutions show considerable ampli- 
tudes near the region of the north pole. This fact indi- 
cates that these solutions also exhibit a character of the 
deformation aligned state. The total angular momentum 
J = 18h consists of J n + = 8h and J p - = lOh. In experi- 
ments there are some discussions about the existence of 
K — 8 + band in osmium [PUlfill and tungsten isotopes |T^ 1 
in connection with an existence of the t-band. The in- 
terpretation of these deformation aligned states requires 
more precise analyses theoretically as well as experimen- 
tally. 

IV. SUMMARY AND CONCLUSIONS 

We have carried out the angular momentum and the 
particle number constrained GCM calculation based on 
the three-dimensional cranked Hartree-Fock-Bogoliubov 
approach. Compared to our previous work we have re- 
fined the force parameters carefully so as to reproduce 
the characteristic features of the yrast states of the mass 
region. With our new set of force parameters the neu- 
tron gap parameter A„ almost vanishes whereas the pro- 
ton gap parameter A p still survives with the size of the 
ground state one at J = 18h. In our calculation the band 
crossing occurs at J = 20h. After choosing the force 
parameter set we have performed the 3D-CHFB calcula- 
tion along the prime meridian on the globe of J = 18h 
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sphere. We know from our previous study that a calcula- 
tion along the prime meridian always gives lower energy 
when compared with the calculation performed on the 
whole surface of the globe. Like the previous calcula- 
tion, we have confirmed the existence of a stable TAR 
potential curve this time also. The TAR curve appears 
at rj) = 24° . Furthermore, here we have found a negative 
potential curve near the north pole. We have performed 
the GCM calculation based on this new potential curve 
and discussed the band structure properties of 182 Os. We 
summarize the present work by stressing following three 
points. 

First point is that by including the constrained terms 
in the Hill- Wheeler equation the stability of the solutions 
is much increased and we can recognize in the GCM solu- 
tions the appearance of the wide range of plateaus even 
in such a high angular momentum region as J — I8h. 
This was not the situation in our previous paper . The 
reason is that we can suppress a drift of the mean value 
of the square of angular momentum value by the con- 
strained term in the Hill- Wheeler equation and prevent 
a contamination by the lower angular momentum com- 
ponents present in the cranking wave function. 

Second point is that from the symmetry property of 
the lowest GCM solution we can expect that this nucleus 
should be a tilted-rotor at J = 187i yrast state. The tilt 
angle is about ip — 12°. As mentioned in section 3.3 the 
TAR potential curve will not be a necessary condition for 
an occurrence of the tilted axis rotation. In our approach 
only the GCM solutions properly predict the characteris- 
tic features of exotic collective motions. In the fourth and 
the sixth solutions we have found an s-band character in 
the PAR region. We have found a sign of multi-band 
mixing in this nucleus. If the s ra -band states, such as 
the fourth or the sixth solutions, come into energetically 
lower positions in the higher angular momentum region, 
we can say that the band mixing fl7| will become a key 
concept in understanding the band structure of the nu- 
cleus. 

The last point is that we have found the new 3D- 
CHFB solutions along the prime meridian which have 
been obtained by the tilt-back calculation at J = 18h 
state. From the characteristic features of the solution, 
we can say that a new type of seesaw vibrational mode 
of the proton and the neutron pairing occurs in the back- 
bending region in 182 0s. This mode accompanies the 
large 7-oscillations. Thus a wobbling motion is expected 
to mediate a band crossing between the g-band and the 
s-band in this nucleus. The GCM calculation which in- 
cludes these new solutions to explore the band structure 
of A ~ 180 mass region will be the next attractive sub- 
ject. 
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